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The effect of binary alloy disorder on the ferromagnetic phases of /-electron materials is studied 
within the periodic Anderson model. We find that disorder in the conduction band can drastically 
enhance the Curie temperature Tc due to an increase of the local /-moment. The effect may be 
explained qualitatively and even quantitatively by a simple theoretical ansatz. The emergence of an 
alloy Kondo insulator at non-integer filling is also pointed out. 
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Materials with /-electrons such as the rare earths (e.g., 
cerium) or actinides (e.g., uranium) exhibit a wealth of 
highly unusual thermodynamic, magnetic and transport 
properties 1]. The minimal microscopic model that can 
account for this diverse physical behavior is the peri- 
odic Anderson model (PAM) which describes a band 
of non-interacting electrons hybridizing with localized, 
interacting /-electrons 2]. Depending on the position 
of the /-electron energy relative to the conducting 
band, and on the strength of the hybridization V and 
the local Coulomb interaction [/, the PAM is able to re- 
produce heavy fermion, intermediate valence, and local 
moment behavior. In the local moment regime and for 
large U the PAM reduces to the so called Kondo lattice 
model, which may be employed as an effective model for 
manganites exhibiting colossal magnetoresistance [3], or 
for diluted magnetic semiconductors which show promise 
for applications in spintronics At low enough tem- 
peratures the PAM also describes magnetically ordered 
phases. While antiferromagnetic order is well-known to 
occur close to half-filling [2] , ferromagnetic solutions are 
found far away from half-filling 5]. Indeed, ferromag- 
netism has been observed in various /-compounds [61] . 

Alloys of /-electron materials also display intriguing 
properties. For example, by changing the stoichiomet- 
ric composition of Ce(Pti_a;Ni:r)2Si2 the systems can be 
tuned from the local moment regime at a; = to the 
intermediate valence regime at a; = 1 Alloying in- 
evitably introduces disorder into the system. In general, 
disorder is expected to reduce the tendency towards fer- 
romagnetic long-range order of the /-electrons and thus 
lower the Curie temperature Tc. On the other hand, 
in certain cases disorder is even known to improve the 
stability of ferromagnetism. For example, disorder in 
the conduction electron band caused by the substitu- 
tion of Rh by Co in URhi_2;Coa;Ge leads to a maxi- 
mum in Tc at X « 0.6 8]. Similarly, a maximum in 
Tc is observed in CeGn2$i\2-xGex dX x ~ 1.5 which may 
be attributed to an enhanced exchange interaction be- 
tween the /-electron moments induced by the diffusive 



motion of the Cu electrons Finally, alloying Ce with 
La as in (Cei_a;Laa;)3Bi4Pt3 introduces disorder into the 
/-electron system, which can trigger a transition from a 
Kondo insulator to a dirty metal jKi] . Clearly disorder is 
an important feature of many /-electron alloys and must 
therefore be included in any comprehensive theoretical 
study of such compounds. 

Previous investigations of the disordered PAM fo- 
cussed, in particular, on the effect of nonmagnetic impu- 
rities on the heavy Fermi liquid or the Kondo insulating 
state [ll|, and on the disorder-driven non- Fermi liquid 
behavior in Kondo alloys [12]. Grenzebach et al. [l3| re- 
cently presented a detailed study of transport properties 
of the disordered PAM within the dynamical mean-field 
theory (DMFT) [Tij together with a thorough discussion 
of the development of the field. The effect of disorder in 
the /-electrons on the ferromagnetic phase was investi- 
gated by Meyer who found that the Curie tempera- 
ture is always reduced. 

In this Letter we report results of a detailed study 
of ferromagnetism in the PAM in the presence of alloy 
disorder. In particular, we show that Tc can be sub- 
stantially enhanced by disorder in the conduction elec- 
trons. We also predict Kondo insulator behavior away 
from half-filling at particular values of the alloy concen- 
tration. Quite generally, the Hamiltonian of the PAM in 
the presence of local disorder has the form 



where (Qo-) and f}^ {fia) are creation (annihilation) 
operators of conduction (c) and localized (/) electrons 
with spin cr at a lattice site i. Here the on-site ener- 
gies e{ and are random variables and V is the local 
hybridization between /- and c-electrons. The hopping 
amplitude of the c-electrons is given by tij . The Coulomb 
interaction U acts only between /-electrons on the same 
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site. 

The alloy will be modeled by a bimodal probability 
distribution function, P{yi) = x5{yi — yo + A*^) + (1 — 
^)^{Vi ~ Vo)i where = e^, e{ are independent, ran- 



dom variables with reference values j/o 
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The 



alloy concentration is characterized by the parameter x 
and the splitting between the atomic levels of the alloy 
components by the energy = A"^, A-^, respectively. 
While the concentration x and energy splitting A^ are. 



in general, independent parameters the values x = 0, 1 
correspond to a non-disordered system even if A^ ^ 0. 
Hence 5^ = x{l — x)/S.y is a natural parameter for the 
disorder strength of alloy disorder. 

The PAM with binary alloy disorder is solved within 
DMFT 1J|, which becomes exact in the limit of infi- 
nite dimensions [l^. In DMFT the disordered PAM is 
mapped onto independent impurities, i.e., for each ran- 
dom variable {y^} the action has the form 
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The function rja-n describes an effective dynamical hy- 
bridization of the c-electrons with the bath. It is the 
same for all random variables {yi} and is determined by 
the self-consistency equations to be discussed below. We 
start with the local matrix Green function 



Gl°^(r;{y.}) 



,(3) 



where Tr is the time-ordering operator. Since the Green 
function ([3]) depends on {yi} it is a random function. 
Here we perform arithmetic averaging, i.e., the averaged 
Green function is given by 



Gl°^M= / Y[dy,P{y,)G':'={T;{y,}). 

{y.} 
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In the absence of interactions one then obtains the re- 
sults of the well-known coherent potential approxima- 
tion (CPA) [l3|. Effects of Anderson localization are ne- 
glected in this case but can be incorporated by employing 
the geometric average [l^. The self-consistency requires 
the averaged local Green function ^ to be the same as 
the lattice Green function, i.e.. 
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The local self-energies appear in the k-integrated Dyson 
equation So-n = G^n ~ Ga-n, where Qa-n is the local Green 
function of the non-interacting bath electrons, with 



V iuJn + Ai - £o - Van 



(6) 



Eqs. ([2][6]) form a general, closed set of equations, which 
determine all local, dynamical correlation functions of 
the disordered PAM. 

To understand the effect of alloy disorder on the 
physics described by the PAM it is instructive to inves- 
tigate the case U = first, since alloy disorder affects a 



hybridized two-band system in several interesting ways. 
To this end we consider the disorder to act only on the 
c-electrons or the /-electrons, respectively. In the case 
of c-electron disorder the diagonal elements of the local 
Green function ^ are given by 
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The case of /-electron disorder is obtained by exchanging 
/ ^ c in Eq. ([7]). Large energy splitting A° leads to 
a band splitting of the conduction electrons as in the 
single band model [l^, i.e., each alloy subband contains 
2xNl and 2(1 — x)Nl states, respectively, where Nl is 
the number of lattice sites. At the same time, the /-level 
does not split. Altogether the alloy with hybridized c- 
and /-electrons can therefore be a band insulator even 
for total densities different from integer values (2 or 4) 

We now include the interaction U between the /- 
electrons and investigate its influence on the alloy sub- 
bands. The effective two-orbital impurity problem in 
the presence of disorder, Eqs. (I2][6]), is solved by finite 
temperature determinant quantum Monte-Carlo (QMC) 
simulations. Ferromagnetic instabilities are detected by 
the divergence of the static spin susceptibility and by a 
non- vanishing value of the magnetization [20|. In the 
numerical examples presented below the DOS for the 
non-interacting c-electrons has the model form Nq^s) = 
V4 — e^/27r, where the energy unit is t = 1. In the fol- 
lowing we fix the interaction and the hybridization at 
U = 1.5 and V = 0.6, respectively, and include disorder 
either in the /-electron or c-electron system. 

As shown in Fig. [T] the computed Curie temperature 
for the transition to the ferromagnetic state is a non- 
monotonic function of the alloy concentration x. In par- 
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ticular, the behavior is quite different for disorder acting 
on the /- or the c-electrons. 



/-electron disorder: In agreement with Meyer 15| the 
presence of /-electron disorder always reduces the Curie 
temperature relative to its non-disordered values at a; = 
or 1. For strong enough disorder eventually vanishes, 
e.g., at a; = 0.28 and x = 0.75, respectively, for A-'' = 1.7 
(left panel of Fig. [T|). This is due to the splitting of the 
/-electron band at large A-'' which increases the double 
occupation of the lower alloy subband; this reduces the 
local moment of the /-electrons and thereby Tc- 

e-elcctron disorder: By contrast, c-electron disorder 
leads to a much more subtle dependence of Tc on con- 
centration X. Namely, for increasing energy splitting A'^ 
there are, in general, three different features observed, 
the physical origin of which will be discussed in more 
detail later: (i) at a; = 1, i.e., in the non-disordered 
case, Tf. is reduced; (ii) a minimum develops in at 
X = TT-tot — 1 > 0; (iii) is enhanced over its non- 
disordered values at a; = or 1. Altogether this leads 
to a global maximum in Tc vs. x. While the decrease 
of Tc at a; = 1 is a simple consequence of the reduction 
of the energy difference between the /-level and the c- 
electron band, e'^ — = eg — — A"^, for increasing A'^, 
the latter effects are more subtle. 

To understand the minimum in Tc vs. x we com- 
puted the evolution of the spectral functions N'^{uj) at 
X = 0.3 for A^ increasing from to 6 (Fig. There 
is an opening of a gap at the chemical potential sig- 
nalling a metal-insulator transition in this system. This 
is caused by the splitting of the c-electron band due to 
binary alloy disorder and the correlations between the /- 
electrons. Namely, for energy splittings A'^ much larger 
than the width of the c-electron band the total number 
of available low-energy states is reduced from 4iVi to 
[4 - 2(1 - x)]Nl = 2(1 + x)Nl, whereby the filling ef- 
fectively increases by a factor of 4/[2(l -|- x)], such that 
nl^^ = 2ntot/(l + x), if ntot < 2(1 + x). For the filling 
ntot = 1.3 studied in the Figs. [T] and [2] the concentra- 
tion X = 0.3 is a special case since then n^^^ = 2. The 
system is then effectively at half-filling and behaves as a 
Kondo insulator at large U , A'^, and low temperatures. 
In particular, itinerant ferromagnetism is unfavorable in 
this case, i.e., Tc = in the vicinity of x = 0.3 at A"^ = 2, 
cf. Fig. [TJ The metal to Kondo insulator transition at 
non-integer filling in the PAM predicted here is a counter- 
part of the Mott-Hubbard metal to insulator transition 
at non-integral fillings in the one-band Hubbard model 
found in 19|, [21 1 . 

We now turn to the maximum in Tc vs. x. It 
can be understood within the following model based 
on an ansatz for the Curie temperature, Tc{U,V, ^) = 
T^{U, V, ^J)F''{^l - e°c)Ff{n - e°j), which implies that the 
formation of local /-electron moments (F-^) is assumed to 
be independent from the c-electron mediated ordering of 
those moments {F'^). In fact, for the RKKY model this 
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FIG. 1: Curie temperature Tc as a function of alloy concen- 
tration X and energy splitting A-'^ (left column) and A'^ (right 
column) for ntot = 1.3 and — e'} = 3.25. Strong c-electron 
disorder enhances Tc compared to its values at a; = or 1. 
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FIG. 2: Spectral function of c-electrons for different A'^ at a; = 
0.3 (other parameters as in Fig.[T} obtained within QMC and 
maximal entropy at T = 1/60. By increasing A'^ a pseudogap 
opens which becomes a real gap for T — > 0. 



ansatz can be microscopically justified within a static 
mean-field theory [22]. The two functions F'^, F^ are 
determined by Tc calculated within DMFT for the non- 
disorder case at fixed M — £§ o'' ~ j respectively; they 
arc shown in Fig.[3l^a) and[3ljb) for one set of parameters. 
The prefactor T^ is determined by the requirement that 
the dimensionless functions F^ and F'^ be equal to one at 
their maxima. We note that T-^(/i — Eg) has a maximum 
when the /-level is half-filled (/a = Eq + U/2), i.e., when 
the local moment is maximal. 

The Curie temperature in the presence of c-electron 
disorder can now be estimated by averaging over the c- 
electron part, T'^, giving rise to the disorder-dependent 
function J^'=(x,/i-eg) = [xF'=(/a-eg-hA'=)-f (l-x)F^(/i- 
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FIG. 3: (a) F^ip - Eq), (b) - e^) appearing in the 

ansatz for Tc (see text) calculated for A'^ — 0. (c) J-'^{fi — Eq) 



and F-^ {fi — Eq ) for A'^ = 2.0; other parameters as in Fig. [T] 
(d) Comparison of Tc obtained from the ansatz and within 
DMFT. 



Eq)]. The linear dependence on the alloy concentration 
can again be justified microscopically within a static 
mean-field theory for the RKKY model, where Tc de- 
pends linearly on the DOS at the chemical potential [2^ . 
Tc is now determined for each concentration x. We cal- 
culate /Lt, which is an implicit function of x, in the non- 
hybridized limit {V = 0) within a rigid band approxi- 



23[. The dependence of the resulting functions 
Eq) on X are shown in Fig.[3jc) 



matron 

J^^(x,/x-eg) and Ff{fi 



for A"^ ~ 2.0. In general {p,~€q) has a global maximum 
at those values of x for which the /-level is half-filled [see 
Fig. [3l[c)]. By contrast, !F'^{x,ii — eg) is characterized by 
a wide minimum, related to the formation of the pseudo- 
gap in the interacting DOS seen in Fig. [51 This minimum 
reaches zero, i.e., J-'^{x,^ — eg) = 0, for a finite range 
of X values as shown in Fig. [3]Jc). The resulting Tc{x) 
obtained by the product of these two functions agrees 
remarkably well with the numerical result obtained by 
DMFT as shown in Fig. [Sid). 

In conclusion, the interplay between the disorder in- 
duced splitting of the conduction band and many-body 
correlations among the /-electrons can lead to a remark- 
able enhancement of the Curie temperature in the peri- 
odic Anderson model. There are two competing effects 
determining as the alloy concentration x is decreased 
from X — 1: (i) a rise due to an increase of the local 
moment, and (ii) a decrease due to the opening of a gap 
in the alloy Kondo insulator at non-integer filling. Al- 
together this leads to a global maximum in Tc vs. x. 
Therefore experimental investigations of /-electron ma- 
terials with alloy disorder in the conducting band are 
expected to be particularly rewarding. 
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